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A recent study by R. L. Liboff of electromagnetic wave propagation in plasmas with a specific 
type of two-dimensional applied magnetic field is continued in this paper. The plasma is assumed to 
be of uniform density in general, and a "spoke-wheel" magnetic field is considered which varies as the 
inverse radius. Perturbation series are obtained for the first Fourier component of the electric field 
for several limiting cases. Both the separate cases of a radial and a circular magnetic field are con- 
sidered in the present work, the former in greater detail, and collisions are included. By suitable 
nondimensionalization, two parameters arise naturally in the present analysis, thus affording a more 
general solution which should be particularly useful in subsequent applications. 

1. Introduction 

In the present paper we shall extend the investigation initiated by Liboff [1964] of the cold 
plasma waves that may propagate in a specific type of two-dimensional magnetic field. A back- 
ground discussion of this problem is given by Liboff. Both the separate cases of a radial and a 
circular magnetic field are considered in the present extension, the former in greater detail, and 
collisions are included. The four physical variables of the problem (electron density, collision 
frequency, magnetic field strength, and geometry) are included in two parameters which arise in 
the present analysis, in which perturbation series are generated for the electric field. The limiting 
cases are suggested naturally on this basis, and quadratures are obtained for the first-order terms 
in the perturbation parameter. The general form of the parameters results in a systematic solu- 
tion which should be useful in subsequent applications. 

2. Analysis 

2.1. Equations 

A collision term has been added to the matrix A in Liboffs formulation, and the governing 
equations were nondimensionalized in a convenient form to facilitate perturbation solutions (refer- 
ence can be made to the detailed presentation [Gold, 1964] for specific intermediate equations and 
results). Assuming solutions of the form 

E a (p, 6)= J EaeWe*** (a = p, 0, z), 

0=0 (1) 

we obtain the following generalized governing equations for the electric field E: 

(± ipip 2 )[pE' etJ + E 9li - (± ifi)Er&] = («or ) 2 {[l ~A(W- Y^Erg + AY r Y e E ee - iAUYeE*} , (2a) 
E'en + -E' m -- 2 Eee - (± WW* + (± #/p*)&* 

= - (nor nAY r YeE rl3 + [1 - A(U* - Y$\E» e + iAUY r E iP } , (2b) 

K& + - E z0-^ E *0 = - (woWA UY e E r0 - iA UY,E m + Q.-A U^E^], (2c ) 



*Present position: Manager, Aerospace Physics Department, Space Systems Division. Hughes Aircraft Co.. Kl Segundo, Calif. 
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where the radial dimension has been nondimensionalized by some reference length r , i.e., p = r/ro, 
and the prime denotes a derivative with respect to p. The following notations are employed: 
A=X/U(U 2 -Y 2 ), Y 2 = Y 2 r +Y 2 e , U=l-i<0 c la>, X = Ne 2 le m(o\ and Y = eB /mw, where w is the 
angular wave frequency, a) c the collision frequency, no = 27rlko the free space wave propagation 
constant, Ao the free space wavelength, eo the free space dielectric constant, and N, e, and m are 
the number density, charge, and mass of the electron, respectively. We shall consider only the 
zeroth mode in the present analysis, and delete the subscript /3 — for convenience. The govern- 
ing equations for the two explicit cases of interest; namely, case 1: Y$ = 0, and case 2: F r = 0, 
follow immediately from equations (2). 

2.2. Solution of Case 1 

The E r field is completely uncoupled from Ee, z in this case and is nontrivial only for co c = 0, 
X=l, i.e., when the signal frequency is equal to the plasma frequency. The remaining two coupled 
equations define Ee and E z . Introducing Y=a/p, these become, 

_ M P { P Ee-iNE z ) Mp( P E z + iNE») ^ 

Llte ^Z^i ' LoEz- pTZTw ' (3a, b) 



where, 



dp 2 p dp 



{n r y 

P. 



(4) 
M = (nor ) 2 XIU, N=a/U. 

The four nondimensional constants of the problem, namely, X, w c /cd, r /Ao, and a combine to give 
two parameters M and TV which must be ordered in the present expansion solution of (3). 

Case la: \N\ > I, i.e., \a\> \U\ = [l + ^/w) 2 ] 1 / 2 . 

Let us assume that the plasma region is annular, n ^ r ^ r 2 , and choose the reference length 
ro = r 2 so that ri/r 2 ^ p ^ 1. The function (p 2 — TV 2 ) -1 can then be expanded in inverse powers of 
TV, and solutions for the zero mode components obtained in the form 

E(p)=jE<»)(p)iV-». 

n=0 

It is now necessary to order M relative to TV systematically. 
(i) \M\ > 0(|7V 2 |). 

The solution for E (n) degenerates into a null series, 
(ii) |M| = 0(|/V|). 

Since \MjN\ is of order one, the zeroth order equations remain coupled. In particular, to 

(XTV- 1 ): 

UEf - ^ £W) = 0, loEja + ^ Ef = 0; (5) 

L ^e N ^z TV * ' z N e TV z ' ^ 

A single fourth-order linear equation can be obtained from (5) and, for example, a power series 
solution derived. Alternatively, power series solutions can be obtained directly from the system 
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(5) and (6). Since the procedure is substantially different than the method of solution used in each 
of the remaining cases, we will not carry out the details at this time. 

(iii) |M| = 0(1). 

The zeroth and first-order equations are 

£,£<<»(£) = 0, B„Ef\{) = 0; (7) 

where the change of independent variable, ^ = «o"oP, introduces the Bessel operator of order n 
and argument £, 

D <P . 1 d ,/. n\ 

B -de + ~^ + V~W (9) 

Hence, 

£»> = cMi) + CiYtf) m Ptf), (10a) 

Eg»=cJ,(£) + C4F,(£) - Pi(f). (10b) 

The solution of (8) and (10) may be obtained using a procedure similar to the one used by Liboff 
[1964]. The particular solutions for E (1) are derived by means of the substitutions, E^df; ) — € z Po, 
E\}'(£;) = €f)P\, where Po = C:Jo+C4Y,> and P\ = cJ\+c-iY\. Thus, 

-M ftff ., . f d£ , 

e * = w^fM d t +c: >lm +CK ' (lla) 

me c di , 
e9= i(w +C7 JiPl +C8 - <llb) 



The solutions for the corresponding homogeneous equations can be absorbed into the zeroth order 
solutions, and one obtains, to order TV -1 , 

E z = P -hN- 1 Po€z, Ee = Pi+N-ip l ee, (12) 

where e is given by (11). Note that in order to satisfy boundary conditions in general two arbi- 
trary constants are required in both the zeroth and first-order solutions; the constants for the latter 
are included in the expressions for e. 

The extension of this solution to include \M\ ^ 0(|7V -1 |) can be readily deduced from the pre- 
ceding results. 

Case lb: \N\ < 1, i.e., \a\ < \U\ = [l + Ww) 2 ] 1 / 2 . 

Let ro=ri in this case so that 1 ^ p ^ r2/n. The function (p 2 — /V 2 ) -1 can now be expanded 
in ascending powers of AT, and solutions for the zero mode components obtained in the form 

E(p)= J E<'<>(p)yV\ 



It is necessary, once again, to order M relative to N. 
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(i) \M\&0(\N- l \). 

The solution for E (rt) degenerates into a null series. 

(ii) |M| = 0(1). 

B 1 E^K0 = 0, BoE?\Q = 0; (13) 

where, a = (l—X/U) il2 is the index of refraction with no magnetic field, and the Bessel operator 
involves the argument £, = noroap. Proceeding as before, one obtains, to order N: 

EM) = Po + NPoe,, E e (Q = P t + NP.ee, (15) 

-iM 



2n r a 



ln£ + c 5 






- iM (H^+ (<% + 



(hi) |M|=0(|JV|). 



£,£<<»(£) = 0, B Ef\i) = 0; (17) 

W-Jjfe*. *W'j£j*. (18) 

The particular solutions of (18) can now be assumed to be of the form, E i z \ ) d = e z ,dE i ® ) d , and one 
obtains: 

E(£) = E«»(1+A4 (19) 

where £»> = P«(©,£i > = Pitf), 

and 

It should be pointed out that a variable electron density could be included in the present 
quadratures, for example (or in the preceding solutions by similar manipulations), in the following 
manner. Let x(£)= X(£)l%m be the normalized radially varying electron density, where X m is the 
maximum value of X in the plasma (hence, x ^ !)• A new constant parameter Mi can then be 
defined in place of M, namely, M = M\X-> where M\ = {n {) r { )) 2 X nt /U. The above solution for e z would 
then become: 

The extension of case lb to include \M\ ^ 0( | TV 2 1 ) can be readily deduced from the preceding 
results. 
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Case lc: |/V|=0(1), i.e., \a\ = 0(|i/|) = 0( Vl + (ov/o>) 2 ). 

Since (p 2 — A^ 2 ) -1 can no longer be expanded, we are left with the single parameter M in (3). 

(i) |Af|>l. 

An expansion in inverse powers of M would result in a degenerate null series for the electric 
field. 

(ii) |M|=0(1). 

No parameter is available for obtaining expansion solutions to (3). It should be remarked, 
however, that only \M\ and \N\ were specified in the discussion of cases la and lb. Thus the 
quantities \a/U\ and \(noro) 2 X/U\ or \{noro) 2 Xla\ were specified. However, the actual magnitudes 
of n ro, a, X, and \U\ were restricted only insofar as they affect these ratios. By further speci- 
fying these parameters, the solutions in cases la and lb can be simplified, say, in evaluating the 
quadratures. Limited solutions can also be obtained in the present case on the same basis. For 
example, let us now nondimensionalize the radial dimension by the free space wave propagation 
constant so that x = n {) r. In place of (3), one obtains the following system of equations: 

D p Mx(xE e - iNE z ) _ _ Mx(xE z ± iNE ) 

B ^= x *-w ' BoE * = * 2 -/v 2 ' {22) 

where Y= a/x, 

and M=XIU, N=a/U. (23) 

Hence, by restricting the analysis to sufficiently large value of \ (for fixed n and r 2 ), so that 
x << 1, expansion solutions in powers of x could be generated. Similarly, for sufficiently small 
value of Xo, x » 1, and expansion solutions in inverse powers of x could be obtained. 

(iii) \M\ < 1, i.e., (n {) r {) ) 2 X < [1 + (<o c /a>) 2 ] 1/2 . 
Assuming solutions in the form 

E(p)= J E<«>(p)M w , (24) 



eq (3) reduce to 



2 n-W" } = -^5 2 (pE^-iNE^)M", 



(25) 

to X Ei " )M " \ = -7^p 2 (pE^ + iNE^M", 

which can be solved, provided that n/r 2 ^ p ^ 1, i.e., by selecting r = r 2 for the reference length. 
General recursion relations can be obtained in this case (/i = 0, 1, 2, 3, . . .): 

(26a) 

^ l£j e „2 AT2 ' °^ ^2_/U2 (ZOO, C) 



£,£«» = 0, 


LoE[» = 0, 


p(pE^-iNE^) 
p 2 -N 2 


pipE^ + iNE"') 



To order M, therefore, 



E(£) = E«»(1 + Me), (27) 
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where 



and 



Ef = fttf), Ef = PM), i = n r oP , 



(n r ) 



1 f dj f eUtPo + inproNPi) f d£_ 

<ro) 2 J £P% J e - (noroN) 2 ** * J ^ C6 ' 

A variable electron density could be included in the manner outlined at the end of the discussion 
of case lb (hi). 

2.3. Solution of Case 2 

In this case E$ is completely uncoupled from E r and E z and is independent of the applied mag- 
netic field Be. Indeed, the following exact solution is immediately apparent: 

& = cJ,(O + c*r,tt)-/»i(0, (29) 

where t, = n,(,r(,p{\—XjU) xl2 . Since the reduced eq (2a) is algebraic, the equations for E r and E z 

can be uncoupled. Assuming Y=—, one obtains 

P 

(iXlaU)Np Mp* 

Er ~ p 2_^2 h ^ Lo/s z (p)- p2 _^ 2 £;, (30a, b) 

M = {norofXIU, N = ^ a=(\-^j' 2 - (3 1) 

The problem reduces to the solution of (30b) for E z in terms of the two explicit parameters, 
M and N. E r can then be obtained immediately from (30a). An approach similar to the one em- 
ployed in case 1 can be adopted in order to solve (30b), although this is not the only possible pro- 
cedure. We will limit our analysis of this problem to the following brief consideration. 

The right-hand side of (30b) can be neglected if \Mp 2 /p 2 — N 2 \ < 1. This condition is satisfied 
if \M\ < 1, ro is selected such that p ^ 1 and the possibility of TV being real with N ~ p is excluded. 
In this case (30b) becomes: 



(norone-nyypy 

and a series solution in powers of Mi can be obtained. A variable electron density has been in- 
cluded as outlined at the end of the discussion of case lb (iii). It follows that: 

EM) = Ei 0, i 1 + Mic), £<»> = Ptf ), (32) 

■d£ f exPl ,,, f di 



(n r y 



[dj r Fxn CM 

){Pl) e-(n rofiY d S + C:i ) tP% 



Cl. 



The corresponding solution for E r is given by (30a). 
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3. Conclusions 

The introductory analysis has been completed on a class of electromagnetic waves which 
will propagate through a cold plasma in which is imbedded a specific type of cylindrical magnetic 
field. The plasma is assumed to be of uniform density and collision frequency, and a so-called 
"spoke-wheel" magnetic field is applied which is both anisotropic and inhomogeneous (varying 
as the inverse radius). Perturbation series were obtained for the first Fourier component of the 
electric field by systematically ordering the two parameters of the problem. The results are given 
in quadrature form and exhibit the correct limiting character in the lowest order expressions. 
Similar studies of guided electromagnetic waves through cylindrical plasmas with imbedded B z 
fields (constant) have been made [Liboff, 1964]. The electric field components are uncoupled 
in these cases, whereas in the present analysis a cylindrical magnetic field is applied and the 
electric field components are severely coupled. 

It remains to fit these solutions to specific geometries. The most physically relevant of these 
is the class of annular geometries referred to in the text. The formulation and solution of a mean- 
ingful boundary value problem based on the present results would enable one to evaluate quan- 
titatively the effect of a variable magnetic field on propagation as a function of the several physical 
parameters of the problem. 
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